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ABSTRACT 

Context. The credibility of an eclipse timing variation (ETV) diagram analysis is investigated for various manifestations of the mass 
transfer and gravitational radiation processes in binary systems. The monotonicity of the period variations and the morphology of the 
respective ETV diagrams are thoroughly explored in both the direct impact and the accretion disk mode of mass transfer, accompanied 
by different types of mass and angular momentum losses (through a hot-spot emission from the gainer and via the L2/L3 points). 
Aims. Our primary objective concerns the traceability of each physical mechanism by means of an ETV diagram analysis. Also, 
possible critical mass ratio values are sought for those transfer modes that involve orbital angular momentum losses strong enough to 
dictate the secular period changes even when highly competitive mechanisms with the opposite direction act simultaneously. 
Methods. The J — P relation that governs the orbital evolution of a binary system is set to provide the exact solution for the period 
and the function expected to represent the subsequent eclipse timing variations. The angular momentum transport is parameterized 
through appropriate empirical relations, which are inferred from semi-analytical ballistic models. Then, we numerically determine 
the minimum temporal range over which a particular mechanism is rendered measurable, as well as the critical mass ratio values that 
signify monotonicity inversion in the period modulations. 

Results. Mass transfer rates comparable to or greater than 10“* Mq yr“* are measurable for typical noise levels of the ETV diagrams, 
regardless of whether the process is conservative. However, the presence of a transient disk around the more massive component 
defines a critical mass ratio (q^ ~ 0.83) above which the period turns out to decrease when still in the conservative regime, rendering 
the measurability of the anticipated variations a much more complicated task. The effects of gravitational radiation proved to be rather 
undetectable, except for systems with physical characteristics that only refer to cataclysmic variables. 

Conclusions. The monotonicity of the period variations and the curvature of the respective ETV diagrams depend strongly on the 
accretion mode and the degree of conservatism of the transfer process. Unlike the hot-spot effects, the Lagrangian points L2 and L3 
support very efficient routes of strong angular momentum loss. It is further shown that escape of mass via the L3 point - when the 
donor is the less massive component - safely provides critical mass ratios above which the period is expected to decrease, no matter 
how intense the process is. 

Key words, binaries: close - accretion, accretion disks - gravitational waves - methods: miscellaneous 


1. Introduction 

It is well known that the information taken directly from the ob¬ 
served orbital evolution of eclipsing binary stars is valuable for 
studying crucial physical mechanisms related to the stellar struc¬ 
ture and evolution. In this direction, the analysis of ETV (eclipse 
timing variations) or, equivalently, of O-C (observed minus cal¬ 
culated) diagrams plays an increasingly important role as a tool 
with which the orbital period variations of eclipsing pairs are dis¬ 
closed (e.g., Hilditch rZOOll Sterken l2005l) . In the framework of 
further mathematical treatment, these variations can be used to 
probe evolutionary processes and estimate the parameters con- 
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nected to the underlying physical mechanisms. In this context, 
the current series of papers aims to further evaluate the efficiency 
of ETV diagrams and to specify their limits as a diagnostic tool 
for long-term observed period variations. 

The incorporated mathematical approach has already been 
presented in detail in an earlier paper (Nanouris et al. l201 II here¬ 
after NKARLl 1). In NKARLl 1, the mathematical procedure is 
implemented for wind-driven mass loss and magnetic braking 
(mainly to estimate the minimum temporal intervals required for 
detecting their signal in ETV diagrams). The morphology and 
the general properties of the ETV diagrams under the separate 
or joint influence of these specific mechanisms were also exam¬ 
ined mostly in late-type detached binaries. 
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In this study we continue and further elaborate on the re¬ 
search by focusing on the effects the following three physical 
mechanisms have on the ETV diagram modulation: (a) conser¬ 
vative mass transfer (MT) through the Lagrangian LI point, (b) 
non-conservative MT accompanied by mass and angular mo¬ 
mentum loss (either as a hot-spot emission from the gainer 
or through the L2 and L3 point), and (c) angular momentum 
loss (AML) due to gravitational radiation. A non-conservative 
MT era will also be referred as “liberal”, as first introduced by 
Warner ( 119781 1 and widely adopted in many subsequent studies 
(see, e.g., van Rensbergen et al. 120081 and references therein). 
Linally, as case studies of particular interest, we examine the ef¬ 
fects that the combined action of the adopted physical mecha¬ 
nisms has on the ETV diagram structure, searching for possible 
critical mass ratio values above or below which the orbital AML 
dictates the secular period changes. 


2. Mathematical procedure 

In this paragraph we briefly outline the basic equations and the 
general methodology and notation developed in NKARLll. In 
particular. Mi, Ri and M 2 , R 2 are the masses and the radii of the 
primary and the secondary component, respectively. Their mass 
ratio q - M 2 /Mi is imposed henceforth to take values below (or 
equal to) unity only, i.e., 0 < q < 1. 

The J - P relation of Kalimeris & Rovithis-Livaniou (120061 1 
is first adopted in order to associate the variations in the state 
of the main physical parameters describing a binary system with 
its orbital period. The relation generalizes the form of Tout & 
Hall (II 99 II 1 by allowing mass losses (ML) from both compo¬ 
nents, also taking the influence of their spin angular momenta to 
the period variations into consideration (i.e., the spin-orbit cou¬ 
pling). 

Throughout our search, we assume for simplicity circular or¬ 
bits (e = 0) with no developing departures from circularity at all 
(e = 0), favoring mostly short- and intermediate-period binaries 
as tidal evolution theory suggests (e.g., Koch & Hrivnak ri981l 
Giuricin et al. I1984I I. Lurthermore, wind-driven ML (including 
magnetic braking for late-type components) are omitted from the 
present analysis for convenience (i.e., M^j - 0, i = 1,2), aim¬ 
ing to constrict our approach to the exchange processes that take 
place in a close pair, reducing at the same time the free parame¬ 
ters of the physical problem. 

According to these assumptions and simplifications, the gen¬ 
eralized J - P relation of Kalimeris & Rovithis-Livaniou (120061 1 
is reduced to 


P 3(l-q),.^ 3q + 2 . ^ 

— =-M 2 - m - 3-3-1- 3-. 

P 3^2 iff 1 + 3^2 J orb 7orb 7orb 


( 1 ) 


On the left-hand side (hereafter l.h.s.) of Eq. ([T]), P - P(t) 
is the orbital period and P = P(t) its rate of change, while on 
the right-hand side (hereafter r.h.s.), the first term refers to the 
MT through the LI point, with M 2 > 0 to declare the rate at 
which mass is accreted by the gainer, supposing that the donor 
is the more massive star. ML that accompanies non-conservative 
manifestations of the MT process is also included with m < 0 to 
express the loss rate at which mass is rejected by the secondary 
component and, subsequently, escapes from the system. 

When the donor is the less massive star, the approach may be 
easily inverted in the opposite direction. Designating as Mi >0 
the rate at which mass is accreted by the gainer and m < 0 the 


loss rate at which mass is rejected by the primary component, 
Eq. (IB is thus modified to the following alternative form: 


P 

P 


3(g-l) 

M2 


■Ml - 


3q ^ + 2 , Ji 72 ^ J 

- m — 3-3-[- 3-. 

Al\ -(- AI 2 7orb 7oi-b 7orb 


( 2 ) 


The third and the fourth terms in Eqs. ([T]) and (|2]) induct 
the spin angular momentum contribution of each individual star 
by means of the consequential - advected - angular momentum 
transport, following the formalism of Gokhale et al. (120071) : i.e., 
7i = jiMi and ^2 = j 2 M 2 - More particularly, in Eq. ([1]), ji and 
j 2 are the specific angular momenta of the matter leaving the 
donor and the matter arriving at the accretor at rates M\ < 0 
and M 2 > 0, respectively, whereas the opposite framework with 
Ml > 0 and M 2 < 0 holds in Eq. dB. 

Although the donor is not considered synchronized with the 
orbital motion in Gokhale et al. ( 120071 ). here it is assumed to be 
tidally locked, following the approach of many earlier studies 
(e.g., Nelemans et al. 120011 Marsh et al. 120041) . Even for the 
gainer, tidal torques are disregarded from our approach owing 
to the wide range of the synchronization timescales, dictated by 
the friction mechanism and the intrinsic properties of the stellar 
outer envelopes (Z ahnlT977l Zahn fT989l Zahn & Bouchet fT989l 
Tassoul & Tassoul fl992al b'). 

In Algols, for instance, tides seem incapable of spinning 
down the gainer during the transfer process (Dervi§oglu et al. 
120101 Deschamps et al. 120131 ). considering that radiative damp¬ 
ing is the driving tidal dissipation mechanism for stars with a 
radiative envelope (Zahn (19771 ). By using the appropriate tidal 
torque and apsidal motion constants of Zahn (119751 see Table 
1 of his work) for the gainer of lO UMa (M\ =2.1 Mq, e.g., 
Soydugan et al. l2013t . for example, a long synchronization time 
of 9.69 X 10^ yr arises. In stars possessing a convective enve¬ 
lope, on the other hand, turbulent viscosity becomes the dom¬ 
inant tidal-braking process, a much more effective mechanism 
than the radiative dissipation (Zahn [T977l ). Considering the low- 
mass accretor of the system RY Aqr (Mi = 1.27 Mq, e.g., Helt 
119871) . which has fractional radius (with respect to the orbital) 
similar to that of the gainer in lO UMa (~ 0.17), a timescale 
that is four orders of magnitude shorter, 5.23 x 10^ yr, is esti¬ 
mated. Tides may also play an important role in the spin rates 
of double white dwarfs and cataclysmic variables. Although the 
viscosity of degenerate matter suggests timescales on the order 
of 10 '^ yr, turbulent dissipation and the excitation of the non- 
radial modes are able to shorten the synchronization time in val¬ 
ues even lower than 500 yr (see Marsh et al. 120041 and references 
therein). As a result, neglecting tidal coupling might affect our 
approach mostly when turbulent convection is the principal fric¬ 
tion process, expected in late-type stars and, possibly, in degen¬ 
erate components. 

In both Eqs. ([TJ and (O, the fifth r.h.s. term represents the 
angular momentum lost by the system (at a rate j), while Jorh 
accounts for the orbital angular momentum: 


7orb — (3) 

In Eq. (O, Mred = Mi M 2 /(Mi - 1 - M 2 ) is the reduced mass of 
the system, Aorb = G'^^(Mi -h M 2 )^^^P^^^I{2ji)^I^ is the orbital 
radius, Gkep = 2nlP is the Keplerian angular velocity, and G the 
gravitational constant. 

The differential Eqs. ([TJ and (|2j can be solved analytically 
for the orbital period function P{t) when the action of each 
individual physical process (or any combination) of our inter¬ 
est is given or even simulated. In consequence, if e denotes a 
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real-valued - continuous - measure of the orbital cycle (which 
also includes integer values when primary light minima are con¬ 
cerned, e.g., Gimenez & Garcia-Pelavo 119831 Borkovits et al. 
120031 1. the orbital period as a function of time, P(t), is converted 
into the orbital cycle domain, /’(e), through the t - f(e) transfor¬ 
mation that emerges by solving for t the following equation; 

df 

- = P{t). (4) 

de 

In this way, the function A/(e) that practically represents the 
synthetic (theoretically anticipated) ETV diagram is eventually 
determined, supposing that the calculated times of minima are 
supplied by a linear ephemeris with initial minimum time Tq 
and reference period P^. Indeed, as the Pit) function is given by 
Eqs. ([T]) or (|2|i along with its cycle transformed counterpart P(e) 
(through Eq. (IH), then AT(e) can be specihed by solving the 
following differential equation (e.g., NKARLl 1); 


dA7’(e) 

de 


P(e) - P.- 


(5) 


All involved ordinary first-order differential equations (i.e., 
Eqs. ©/©I and (|5])) are solved under the initial condi¬ 
tions (e,t, P(e), AT{e)) — (0, Tq, Pe,0) that eventually provide 
the A7’(e) function. If, additionally, e denotes the noise level of 
an ETV diagram, the critical number of epochs emin, needed for 
the detection of a signal induced by a certain mechanism, is then 
made available by solving the equation. 


A7’(e) = e. 


( 6 ) 


which, alternatively, provides the number of epochs such that the 
signal-to-noise ratio equals one. 

Similarly to NKARLl 1, the value of e = 0.001 d is con¬ 
sidered here as a typical noise level associated with photoelec¬ 
tric times of minima, e = 0.03 d is adopted as a typical value 
associated with visual or patrol plates times of minima, while 
s - 0.01 d serves as a pilot statistical noise level (i.e., for cases 
where the presence of migrating spots cannot be excluded, see 
also Kalimeris et al. I2002I I. 


3. Examination of physicai mechanisms 

The implementation of the preceding methodology follows in 
four sections. In the hrst we examine the effects that MT through 
the LI point (including those induced by a hot-spot emission) 
has on ETV diagrams, while the second section considers the 
effects of ML through the outer Lagrangian points L2 and L3. 
In the third section, we design non-conservative MT schemes 
that unveil the ETV diagram curvature based on the mass ratio 
of a system for a certain degree of conservation. Einally, in a 
fourth, the impact of AML via gravitational radiation in the O- 
C modulation, is examined. 


The inflowing matter may either hit the gainer directly (here¬ 
after direct-impact mode), or it can orbit around the gainer, pro¬ 
ducing internal particle collisions due to intersecting trajectories. 
A gaseous ring is built in this way and eventually turns into an 
accretion disk owing to the effects of the gaseous material bulk 
viscosity (hereafter disk formation mode, e.g., Kruszewski l 1 9671 
Kffz ll97n Lubow & Shu ll975l Lubow & Shu ll976l Papaloizou 
& Pringle [19771 Lin & Papaloizou [T979J Shu & Lubow [T981l l. 
The latter occurs when the detached component has smaller ra¬ 
dius than a critical value, Rmin = ^min^orb, with rmin to be given 
as a function of the mass ratio of the system; 

• = 0.04838 -b 0.03601 log^ -b 0.01229(logg)^ (7) 

when the donor is the more massive component, and 

• = 0.04930 - 0.03387 logq + 0.05915(logqf, (8) 

when the donor is the less massive component. 

This approximation has emerged from the ballistic calcula¬ 
tions of Lubow & Shu (119751 see Table 2 of their work) and is 
valid for 0.07 < q' < 1 with accuracy close to 1% when the donor 
is the more massive component and better than 0.5% when the 
donor is the less massive component. 


3.1.1. Mass transfer in the absence of an accretion disk 

In the direct-impact MT mode, the flow stream coming from the 
donor hits the gainer and dissipates without creating an accretion 
disk when the radius of the gainer is much larger than the respec¬ 
tive critical value Rmin (e.g., Kffz I197T1 Lubow & Shu I1975l l. 
The stream strikes the gainer almost perpendicular to its surface 
and generates enough turbulence to keep the friction timescale 
short (Biermann and Hall 1 19731 Kaitchuck et al. ll985t . As a re¬ 
sult, tidal interaction seems efficient enough to keep both com¬ 
ponents in a strong coupling regime, canceling the angular mo¬ 
mentum flows within the binary (e.g.. Hut and Paczvnski ll9841 
Verbunt and Rappaport I 19881 1, i.e. 72 ~ --/i- The orbital angular 
momentum then remains invariable, and the period varies as a 
result of the mass redistribution between the two components. 

In the fully conservative framework {m = 7 = 0), the J - P 
relation given by Eqs. ([T]i or © is reduced to an ordinary first- 
order differential equation of the form, 

P = WLlP, (9) 


where, when the donor is the more massive star, wli is given by 


. WLI =-3(1-^)^ <0, (10) 

Ml 

while, when the donor is the less massive star, wli is given by 


wli=-3(^-1)^>0. 


( 11 ) 


3.1. Case I: Mass exchange through the LI point 

A semi-detached configuration is attained as long as one of the 
two components of a binary system reaches the limiting surface 
of the Roche lobe. The MT process, also known as Roche lobe 
overflow (RLOE), is then realized through a critical (i.e., gravita¬ 
tionally neutral) point that coincides with the Lagrangian point 
LI for a fully synchronized system (e.g., Kuiper [79411 Kopal 
119561 Plavec [T958l Kruszewski 1 1 963t . 


Similar to NKARLl 1, the coefficient of proportionality be¬ 
tween the rate at which period changes and the period itself, wli, 
is considered constant since M, terms are small enough to af¬ 
fect the absolute parameters of any of the components during 
the time windows covered even by the longest available ETV di¬ 
agrams. Therefore, solution of Eq. (I9]l becomes trivial and under 
the initial condition (f, P{t)) = (Tq, P^) is written as 

P(t) = _ (12) 
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Then, integration of Eq. (|4li under the initial condition (e, f) = 
(0, To) yields 


Table 1. Minimum time spans at e = 0.01 d in absence of an 
accretion disk (conservative case). 


f-ro = -—In(l-WLiPee). (13) 

wli 

As a result, through Eqs. ([nil and (foi l, the orbital period 
P{e) as a function of the continuous time variable e is 


Pie)^ 


1 - WLlEef’ 


(14) 


Subsequently, integration of Eq. (|3 under the initial condition 
(e, A7’(e)) = (0,0) eventually gives the corresponding solution 
for the AT(e) function: 


q 

M. 

[Me yr 

= 0.5 

1 — To 
‘] [yr] 

M 

q 

= 10 - 

^2 

* Me yr 

r\ 

t-To 

[yr] 

KF^ 

6042 

0.9 

0.19 

0.20 

135 

10-" 

1911 

0.8 

0.18 

0.21 

95.5 

lO-'o 

604 

0.7 

0.18 

0.22 

78.0 

10-5’ 

191 

0.6 

0.17 

0.24 

67.6 

10-** 

60.4 

0.5 

0.17 

0.25 

60.4 

10-’ 

19.1 

0.4 

0.16 

0.28 

55.2 

10^^ 

6.0 

0.3 

0.15 

0.31 

51.1 

10-^ 

1.9 

0.2 

0.13 

0.36 

47.8 

10-^ 

0.6 

0.1 

0.11 

0.45 

45.0 


A^e) = —ln(l-WLiPee)-^’ee- (15) 

wli 

Equation (ITSl l implies that the curvature of the AT{e) func¬ 
tion, and subsequently the anticipated form of the respective 
ETV diagram, is determined by the sign of the constant wli. 
In particular, when the donor is the more massive star (wli < 0, 
see Eq. (ITOl ll. the orbital period decreases (see Eq. (j^l) so the 
AT{e) function is concave for any e (see Eq. (fTTT i). Likewise, 
when the donor is the less massive star (wli > 0), the orbital 
period increases while the AT(e) function is convex for any e. 
Therefore, the foregoing analysis recovers the widely used pro¬ 
cedure for estimating MT rates under conservative conditions 
in semi-detached and contact binaries (e.g., Kwee [19581 Huang 
119631 Kruszewski |1964b1 l. 

Given the AT{e) function, Eq. (j^ is solved by a suitable 
Newton-Raphson numerical scheme (e.g.. Press et al. 119921) in 
order to determine the minimum epoch e^m and the subsequent 
temporal intervals f(emin) - To from Eq. (HI. Our test case con¬ 
cerns a ^ = 0.5 semi-detached binary with a main-sequence sec¬ 
ondary gainer of fractional radius r 2 - RilAorh = 0.17 or a main- 
sequence primary gainer of fractional radius r\ - Ri/Aorb = 
0.25, keeping the secondary as massive as before. Results for 
several MT rates with a pilot noise level of s = 0.01 d (in Eq. (0) 
are given in Table [1] while for e - 0.001 d and 0.03 d, the ex¬ 
pected intervals are depicted in Fig. [T^. Similar computations 
were also made for various q values at M\ - 10“* Mq yr“’ (see 
Fig. [TJ)), always taking into consideration that r 2 , ri » rmin 
should be valid. The temporal intervals are expected to be the 
same regardless of the component from which the mass begins 
to flow (given that |1 - ^| remains constant) and that they do 
not depend on the reference period as long as Eq. (0 holds (see 
NKARLll). 

Because most of the available times of minima span slightly 
more than a century, a close inspection of Fig.[T]shows that under 
a low noise level of e = 0.001 d, the signal of the conservative 
MT process without disk formation can be detected at rates M, 
greater than 10“'°Moyr“*, and mass ratios even close to unity. 
For higher noise levels, such as e = 0.01 d or 0.03 d, the signal 
is rendered measurable only when M, > lO^^Moyr ' and q < 
0.9 (see Table[T]and Fig.[TJ. The aforementioned ascertainments 
are in fairly good agreement with MT rates derived so far by 
traditional approaches for different binary types with at least one 
component filling its Roche lobe, ranging from 10“* Mq yr*' to 
10“^ Moyr ' (e.g., Kaszas et al. 119981 Simon [T9991 Qian & Liu 
120001 Qian 1200^ Qian & Boonrucksar 120021 Qian et al. 120071 
Manzoori 120081 Qian et al. 120081 Zasche et al. 120081 Zhu et al. 
120091 Christopoulou et al. l201 111 . 


3.1.2. Mass transfer in the presence of a transient accretion 
disk 

Still in the direct-impact MT mode, the stream may be deflected 
by the Coriolis force before striking the gainer when its radius 
is marginally above the respective critical value Rmin (e.g., Kfiz 
119711 Lubow & Shu [1975). The stream hits the gainer tangen¬ 
tially to its surface, torquing the star near its limb and producing 
a thin ring that is gradually transformed into a quasi-stationary 
disk structure (Kaitchuck et al. 119851 Peters 119891 Richards & 
Albright [19991 l. The disk will be considered as transient (or low 
viscous) when the bulk viscosity is not efficient enough to re¬ 
turn the proffered angular momentum to the orbit (we adopt the 
notation of Kaitchuck et al. 119851 1 or, equivalently, when it ex¬ 
periences a weak coupling regime with the donor. The advected 
angular momentum progressively accumulates in the disk and 
spins up the gainer, forcing the orbital angular momentum to de¬ 
crease as a result of the conservative status (see, for instance. 
Packet [19811 for an excellent review on this topic). 

In the present work, the specific angular momentum arriving 
at the accretor is parameterized in terms of the circularization ra¬ 
dius (e.g., Biermann & Hall 119731 Flannerv ll975l Warner [T9781 
Hut and Paczvnski 119841 Verbunt & Rappaport 11 9881 1, i.e. the 
equivalent radius of a ring — rrAmb in a circular Keplerian 
orbit around the gainer that has the same specific angular mo¬ 
mentum as the transferred mass. Hence, J 2 = and 

ji = VGiTO are the specific angular momenta accumulated 
by the accretor as secondary or primary, respectively, resulting 
in angular momentum gain at rate J 2 IJorb — yjqiq + l)r^M 2 lM 2 
and j\ /7orb = yl(q+ IM 2 for each corresponding case. 

Here, we provide an empirical relation for the fractional ra¬ 
dius of the ring as a function of the mass ratio q, inferred 
from the results of the semi-analytical approach of Lubow & 
Shu (119751) for sufficiently low initial flow stream velocities at 
the LI point (see Table 2 of their work): 

• 0.08586 -H 0.06214log^ H- 0.02022(log^)^ (16) 

when the donor is the more massive component, and 

• r, = 0.08658 - 0.06925 log^ h- 0.07368(log^)^ (17) 

when the donor is the less massive component. 

This approximation is valid for 0.07 <q< 1, with the av¬ 
erage relative error of r, being estimated close to 1% when the 
donor is the more massive component and less than 0.5% when 
the less massive. In the latter case, we also refer to the even more 
accurate relation of Verbunt & Rappaport ( 119881 1 when a wider 
mass ratio range is needed (i.e., for 0.001 < q <\). 
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Fig. 1. Minimum time spans in a conservative era a) for several mass transfer rates when the mass ratio is fixed at ^ = 0.5 and b) for 
various mass ratios when the mass transfer rate is fixed at Mq yr“'. Computations were performed both in the absence (circles) 
and in the presence of a transient disk around either the less massive (squares) or the more massive component (up triangles), for 
noise levels equal to 0.001 d (left vertical axis) and 0.03 d (right vertical axis). The dashed line in the right panel covers the range 
within which the critical value ~ 0.83 lies. (This figure is available in color in the electronic form.) 


As far as the specific angular momentum leaving the donor 
is concerned, if /?l = ^LAorb is the (equivalent) radius of its 
Roche lobe, for the purposes of the present study we simply 
adopt i\ R^Qkep in Eq. ([T]i and 72 ~ in Eq. (|2]), 

following the suggestion of Gokhale et al. (120071 1 and recalling 
that the donor here is synchronized with the orbital velocity and 
fills its Roche lobe. As a consequence, its inertia offers an extra 
store of angular momentum, leaving the donor at a rate equal to 
7 i/ 7 orb = -iq+ l)rlM2lM2 and /a/^orb = -(q + l)rlMi/M2 for 
each case, respectively, considering that Mi = -M 2 in the con¬ 
servative framework. The fractional Roche lobe radius tl is es¬ 
timated via the empirical law of Eggleton (119831 1 with accuracy 
better than 1%, properly modified so that q < I (e.g., Hilditch 

IMB: 

^ + \n{\ + q^^l^y ^ ^ 

when the donor is the more massive component, and 

0.49^2/3 

^ Q.6q^^^+ \n{\ + q^iy ^ ^ 

when the donor is the less massive component. 

As a consequence, the j - P relation, given by Eqs. dB 
and (l2]l, turns into an ordinary first-order differential equation 
identical to Eq. (|9]l, P - wli.iR, which can be solved for the or¬ 
bital period P{t), by taking into account that, when the donor is 
the more massive star, the constant WLi.r is given by 

• wuj = -3 [1 - ^ -I- ^Jq(qTTj^\- -(q+ l)r£] ^ < 0, (20) 
while, when the donor is the less massive star, WLi.r is given by 

• WLi,r = -3 - 1 -I- y(q+ \)r, -{q+ l)rl\ (21) 


Since the presence of a transient disk drains orbital angu¬ 
lar momentum, when the donor is the more massive component, 
the orbital period still decreases with respect to the no-disk case 
{wLi,r < 0, see Eq. (1201 0. However, the sign of WLi.r is not steady 
when the donor is the less massive star (see Eq. (12110 . so the pe¬ 
riod varies according to the degree the viscosity of the transient 
disk facilitates the orbital AML. High viscosity values favor a 
strong gainer-disk-orbit coupling regime that returns most of the 
transferred angular momentum back to the orbit; i.e., wli,i > 0 
and P > 0, while low viscosity values favor a weak coupling 
regime that is inadequate for restraining the orbital AML; i.e., 
wli,i < 0 and P < 0. The extra deposit of angular momentum 
stemming from the spin contribution of the donor proves to be 
too weak to inhibit the orbital shrinkage. 

In the simplified framework of Eqs. (flSl) . (flTl i. (fl^ . and (fl^ . 
the sign of wu r is determined solely by the mass ratio of the 
system q. In particular, the unique root at q = 0.83 of the term 
q - 1 + y(q -I- l)rr - (q + l)r^ inside the brackets in Eq. (l2T1 l. 
defines a critical value q^^ of the mass ratio q so that 

• when q < qa - 0.83, then P > 0 and AT{e) is convex; 

• when q > qci - 0.83, then P <0 and AT{e) is concave. 

The detection of a critical value q^r for the mass ratio of¬ 
fers strong support to concave modulated ETV diagrams in high- 
mass-ratio, semi-detached binaries when the less massive com¬ 
ponent is the donor, contradicting the predictions of the conser¬ 
vative MT process in the absence of a disk. 

Determination of the minimum time interval required in or¬ 
der to be the signal of the examined process observable in an 
ETV diagram initially concerns aq - 0.5 semi-detached binary 
with a main-sequence gainer of fractional radius r 2 = R 2 /Aorb = 
0.07 or ri = Ri/Aprb = 0.09 as the secondary or primary com¬ 
ponents, respectively. Were computed several MT rates for both 
the pilot noise level s = 0.01 d (see Tables |2] and |3]l and the more 
extreme cases of e = 0.001 d and 0.03 d (see Fig.[T^). The same 
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Table 2. Minimum time spans at e = 0.01 d in the presence of a 
transient disk around the less massive component (conservative 
case). 


q 

M 2 

[Mo yr- 

= 0.5 

t — To 

'] [yr] 

q 

M 2 = 

^min 

lO-** Mq 

r2 

yr ^ 

Ti- 

t-To 

[yr] 


6462 

0.9 

0.05 

0.07 

0.08 

97.2 

10 -" 

2044 

0.8 

0.05 

0.07 

0.08 

84.8 

lO-'o 

646 

0.7 

0.04 

0.07 

0.08 

76.0 

lo-® 

204 

0.6 

0.04 

0.07 

0.07 

69.5 

10 -* 

64.6 

0.5 

0.04 

0.07 

0.07 

64.6 

10 -^ 

20.4 

0.4 

0.04 

0.06 

0.06 

60.8 

10 -'’ 

6.5 

0.3 

0.03 

0.06 

0.06 

57.9 

10 -^ 

2.0 

0.2 

0.03 

0.05 

0.05 

55.8 

10 -“ 

0.6 

0.1 

0.02 

0.04 

0.04 

55.1 


Table 3. Minimum time spans at e = 0.01 d in the presence of a 
transient disk around the more massive component (conservative 
case). 


q 

Ml 

[Mo yr- 

= 0.5 

t — To 
'] [yr] 

q 

Ml = 

^ min 

lO-** Mq 
ri 

yr ^ 

Ti- 

t-To 

[yr] 

Ttf^ 

8709 

0.9 

0.05 

0.07 

0.09 

185 

10 -" 

2754 

0.8 

0.05 

0.07 

0.09 

290 

lO-'o 

871 

0.7 

0.06 

0.08 

0.10 

138 

lo-'* 

275 

0.6 

0.06 

0.08 

0.11 

104 

lo-** 

87.1 

0.5 

0.06 

0.09 

0.11 

87.1 

10 -^ 

27.5 

0.4 

0.07 

0.09 

0.13 

76.9 

10 -'’ 

8.7 

0.3 

0.08 

0.11 

0.14 

70.1 

10-5 

2.8 

0.2 

0.10 

0.12 

0.17 

65.7 

10 -“ 

0.9 

0.1 

0.14 

0.15 

0.23 

64.1 


procedure was repeated for various q values at a steady MT rate 
of 10“^ Mq yr“' (see Fig.[TJi) by evaluating that the gainer’s frac¬ 
tional radius always lies between rmin and quantities that are 
also presented in Tables|2]and|3] 

Although the inferred intervals f(emm) - To are found to be 
slightly different with respect to those in the absence of a disk 
when the donor is the primary component, they are severely 
broadened in the opposing framework, requiring MT rates en¬ 
hanced by at least one order of magnitude (i.e., greater than 10“® 
for £ = 0.001 d and 10“* Mq yr * for e - 0.03 d). Beside this, the 
existence of a critical mass ratio q„ « 0.83, around which the 
period is expected to remain approximately invariable, makes 
detecting such a signal in an ETV diagram even harder. In prac¬ 
tice, systems with mass ratio between 0.8 and 0.9 are expected to 
display no visible O-C variations with any of the adopted noise 
levels. 

3.1.3. Mass transfer in the presence of a permanent 
accretion disk 

In the disk formation MT mode, the stream misses completely 
the gainer when its radius is below the respective critical value 
^min (e.g., Kffz ll97Tl Lubow & Shu ll975l l. As a consequence, 
a rotating ring is formed that spreads out and is progressively 
transmuted into a stationary disk structure (e.g., Papaloizou & 
Pringle 119771 Lin & Papaloizou 119791 . The disk will be con¬ 
sidered as permanent (or high viscous) provided that the fluid’s 
material bulk viscosity is efficient enough to return the proffered 
angular momentum back to the orbit (once more, we adopt the 
notation of Kaitchuck et al. 119851 or, equivalently, when it ex¬ 


periences a strong coupling regime with the donor (e.g.. Hut and 
Paczvnski ll9841 Verbunt and Rapr)ar)ort ll9881 . However, mate¬ 
rial at the inner edge of the disk still has angular momentum that 
acts to spin up the accretor at the expense of the orbital angular 
momentum (e.g.. Packet 119811 Marsh et al. 120041 Dervi§oglu et 
al. 120101 Deschamps et al. 120131 . 

The description of the period evolution of such systems does 
not differ at all from those systems that accommodate a transient 
disk, by simply replacing by r 2 - R 2 Morb or r\ - Ri/Amb 
when the gainer is the less or the more massive star, respectively 
(Marsh et al. 120041 . The minimum temporal intervals f(emin) -Tq 
derived in the presence of a transient disk should be considered 
mostly as an upper limit (see Tables |2] and |3] and Fig. [TJ, while 
the results inferred in the absence of a disk should mostly be 
regarded as a lower limit (see Table [1] and Fig.lTJ, reflecting the 
case where r 2 ^ 1 or ri « 1 when the gainer is the secondary 
or the primary, respectively. 


3.1.4. Mass and angular momentum loss through a hot-spot 
emission 


When a semi-detached binary undergoes a MT phase, a hot spot 
is expected to be created in the area that the inflowing material 
impacts the mass-gaining component (whether an accretion disk 
is present or not, e.g., Kruszewski ll964al Smak I 19711 Warner & 
Nather lT9711 Lubow & Shu ll975l Lubow & Shu ll976T l. As long 
as the radiative energy of such a hot spot, strengthened due to the 
limited accreted area (along with the rotational kinetic energy), 
surmounts the gravitational binding energy, the system is subject 
to a short non-conservative era, generally expected soon after 
the onset of the RLOF phase. It is proved that this condition is 
fulhlled provided that a critical MT rate has been reached (van 
Rensbergen et al. 12008t . As a result, a hot spot can be a key 
factor in impelling a binary to a non-conservative phase. 

Here, the AML process is parameterized in terms of the or¬ 
bital specihc angular momentum by adopting the formalism of 
Rappaport et al. (119821 1: 

7hs “ (^2) 


where mbs < 0 denotes the ML rate during the non-conservative 
era and y'hs the - dimensionless - specific angular momentum 
taken away from the orbit of the accreting star; i.e., Jhs = (1/1 + 
q)^, if this is the secondary, and y'hs = (?/l + q )^7 if this is the 
primary star (de Mink et al. 120071 Rensbergen et al. l201 111 . The 
AML term in Eqs. ([T]i and ^ is then equal to 7hs/7oib = 7hs(? + 
l)mhs/M2. 

Attempting to isolate the ML and AML effects caused by the 
presence of a hot spot on the orbital evolution, we proceed by 
disregarding the MT term (i.e., M; = 0), and the spin angular 
momentum of both components as well (i.e., J 2 ~ -J\)\ in other 
words, we explore the impact that a fully liberal era (in which a 
hot spot emits all the incoming material) has on an ETV diagram 
morphology in the absence of a disk. The differential equation 
that governs the orbital evolution of the system, as derived from 
Eqs. ([Ill and (0, therefore has a similar form to Eq. (0, F = 
WLi.hsF, with the WLi.hs constant being given by the following 
relations: 


• VVLl,hs 


-3g^ -2q + 3 

^ -H 1 Al2 


(23) 


when the donor is the more massive component, and 
3q^ -2q-3 


WLl.bs 


^ -H 1 


^^>0, 

M 2 


(24) 
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Table 4. Minimum time spans at s = 0.01 d when a hot spot, 
located on the less massive component, emits all the incoming 
material (fully liberal case). 


9 = 0.5 

mhs t - To 

[Moyr-‘] [yr] 

%s = 

q 

-10-’*Moyr-' 
n t-To 

[yr] 

-10-*' 

8106 

0.9 

0.19 

92.0 

-10-'* 

2563 

0.8 

0.18 

138 

-io-“> 

811 

0.7 

0.18 

268 

-10-® 

256 

0.6 

0.17 

110 

-10-** 

81.1 

0.5 

0.17 

81.1 

-10-’ 

25.6 

0.4 

0.16 

66.8 

-10-" 

8.1 

0.3 

0.15 

57.8 

-10-" 

2.6 

0.2 

0.13 

51.5 

-lo-® 

0.8 

0.1 

0.11 

46.6 


Table 5. Minimum time spans at s = 0.01 d when a hot spot, 
located on the more massive component, emits all the incoming 
material (fully liberal case). 


9 = 0.5 

mhs t - To 

[Moyr-‘] [yr] 

Whs = 

q 

-10-’*Moyr-' 
r\ t-To 

[yr] 

-10-*' 

5027 

0.9 

0.20 

66.3 

-10-" 

1590 

0.8 

0.21 

60.6 

-io-“> 

503 

0.7 

0.22 

56.4 

-10-® 

159 

0.6 

0.24 

53.0 

-10-** 

50.3 

0.5 

0.25 

50.3 

-10-’ 

15.9 

0.4 

0.28 

48.1 

-10-" 

5.0 

0.3 

0.31 

46.2 

-10-" 

1.6 

0.2 

0.36 

44.8 

-lo-'* 

0.5 

0.1 

0.45 

43.6 


when the donor is the less massive component. 

Equations (|2^ and (l24li make it clear that the prevalence of 
the driving mechanism (the ML process itself against the ML- 
induced torque through which the system loses angular momen¬ 
tum) depends directly on the mass ratio of the system. The sign 
investigation of WLi.hs in Eq- (l2?t shows that for binaries whose 
the donor is the primary star, the orbital period is expected to 
increase when q - > 0.72, otherwise the period decreases. 

In contrast, in binaries whose donor is the secondary star, the or¬ 
bital period is expected to increase regardless of the mass ratio 
q values. The above inference clearly shows that the AML effect 
hardly compensates for the ML process, resulting in an increas¬ 
ing period in most cases. 

Determination of f(emin) - Tq takes place following the stan¬ 
dard procedure for the test cases examined in the direct-impact 
mode without a transient disk. The estimated temporal windows 
are quoted in Tables|4]and|5]for s = 0.01 d, while they are graph¬ 
ically displayed in Lig.l^for s = 0.001 d and 0.03 d. Peering at 
the tabulated values, one notes that an O-C modulation is ren¬ 
dered observable for ML rates that exceed -10“^ Mq yr“', barely 
favoring the case of the secondary star as the donor. Moreover, 
a non-conservative era due to a hot-spot emission from the sec¬ 
ondary is expected to be measurable for almost any mass ra¬ 
tio except for a very narrow range between 0.7 and 0.8 (where 
q^i ~ 0.72 lies). 

By reducing the noise level at s = 0.001 d, ML rates stronger 
than -10“*°Moyr“' are still detectable for a century, while no 
significant changes occur by worsening the noise contamination 
at e = 0.03 d (Lig. |^). However, by keeping the ML rate fixed 
at ffihs = -lO^^Moyr^', the O-C variations are visible for e 


- 0.001 d regardless of the mass ratio value, whereas a noise 
level increase at e = 0.03 d limits the traceability of the examined 
mechanism only for very small mass ratios when the donor is 
expected to be the primary star (Lig.|2j5). 

3.2. Case II: Mass and angular momentum loss through the 
L2/L3 points 

In the course of a violent MT process in a semi-detached binary 
(e.g., in super-Eddington accretion flows), mass and angular mo¬ 
mentum may be lost from the system either through the L2 or 
the L3 point when the donor is the more or the less massive 
component, respectively. These are the paths that matter can es¬ 
cape most easily from the gravitational field of the binary since 
it overcomes the lowest potential barrier, as a result requiring 
the lowest energy of all. The formation process and the observed 
features of a circumbinary envelope due to outflows from these 
two points have been confirmed or/and explored widely in litera¬ 
ture (Kuiper [l94U Shu et al. 119791 Soberman et al. 11 9971 Spruit 
& Taam l20()n Taam & Spruit l20OTl Chen et al. 120061 D’Souza et 
al. [2006l Motl et al. [2007l Sytov et al. l2007t Sytov et al. l2009al b: 
Zhilkin & Bisikalo[2009j Bisikalo r20101 Mennickent et al. 120101 
Mennickent et al. 1201211 . 

Here we adopt the AML parametrization proposed by Shu et 
al. (119791 1. which expresses the angular momentum taken away 
from the system via the L2 and L3 points in terms of the or¬ 
bital specific angular momentum, similarly to the hot spot case. 
If mL 2 < 0, mL 3 < 0 defines the ML rate of the escaping matter 
from the L2 and L3 points with - dimensionless - specific angu¬ 
lar momenta ji, 2 , 7 l 3 , then 7 l 2 and 7 l 3 are estimated through 
Eq. (l22l l by properly substituting Whs and jhs for each corre¬ 
sponding case. 

Intending to provide an efficient empirical relation for yL 2 
as a function of the mass ratio q, we employ the values of the 
semi-analytical ballistic model of Shu et al. (119791) . as tabulated 
in Table 2 of their study. The values account for the terminal spe¬ 
cific angular momentum (at very large distances from the system 
where the flow Anally obtains its terminal kinematical charac¬ 
teristics), also taking the oblique impact angle the outflowing 
stream follows into account. In the absence of analytic calcula¬ 
tions for 7 l 3 , we assume that matter leaves the system at the L3 
point in a direction vertical to the line connecting the compo¬ 
nents (thus providing an upper AML limit); i.e., jui - with 
xl 3 denoting the relative position of L3 in the Roche geometry 
coordinates system: 

• ;l 2 = 1-70166-0.49695 log^-0.41981(log^)^ (25) 

when ML and AML occur through L2, and 

• ;l3 = 1.43988 -h 0.55460log^ h- 0.19061(log^)^ (26) 

when ML and AML occur through L3. 

This approximation is valid for 0.05 < ^ < 1 with accuracy 
better than 0.5% for both jui and however, 7 l 3 represents 
the initial specific angular momentum (i.e., at L3), and a more 
reliable error close to 15% should be considered (as estimated 
by comparing the values of initial and terminal specific angu¬ 
lar momentum for jui, which are listed in Table 2 of Shu et al. 

mm . 

Even after taking all assumptions mentioned in the hot spot 
case into consideration and specifying that the AML terms 
juilJmh and 7L3/7orb in Eqs. ([TJ and (|2]) are equal to Jl 2 (? + 
l^)mL 2 lM 2 and jLiiq + l.)mijilM 2 , respectively, the orbital evolu¬ 
tion of the system is then still described by Eq. (|9]l by setting the 
wl 2 constant in place of wli as 


7 














































N. Nanouris et al.: Efficiency of ETV diagrams as diagnostic tools for long-term period variations 


-10 -9 -8 -7 -6 -5 
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(b) 


Fig. 2. Minimum time spans in a fully liberal era a) for several mass loss rates when the mass ratio is fixed ai. q — 0.5 and b) 
for various mass ratios when the mass loss rate is fixed at -lO^^Moyr”'. Computations were performed considering that the 
transferred material is lost from the system through a hot-spot emission, located either on the less massive (squares) or the more 
massive component (circles), and via the L2 (up triangles) or the L3 point (down triangles), for noise levels equal to 0.001 d (left 
vertical axis) and 0.03 d (right vertical axis). The dashed line in the right panel covers the range within which the critical value 
q'cr ~ 0.72 lies. (This figure is available in color in the electronic form.) 


3;l 2(? + 1)^ - -2q' mL2 „ 

• Wl2 = -;- — < 0, 

q \ A /2 

when ML and AML occur through L2, and 
37l 3(?+l)^-2g'-3 mL3 


Wl3 - 


q + I 


■ — < 0 , 
M 2 


(27) 


(28) 


when ML and AML occur through L3. 

It is noteworthy that the sign of the numerator appeared on 
the r.h.s. of Eqs. (l27l) and (1281) is positive for any mass ratio q 
(since jiji, 7 l 3 > !)■ Unlike the hot spot case, specific angular 
momentum is carried away from the system at a rather distant 
point (with respect to the center of mass), leading to markedly 
different AML rates and warranting period reduction (y'hs < !)■ 
As a result, the outer Lagrangian points L2 and L3 turn out to 
be routes through which the excess of matter, surviving the MT 
process, travels out of the system with prominent AML rates, 
ensuring orbital shrinkage. 

To examine the measurability of the ML mechanism through 
the L2 and the L3 points, the standard methodology followed 
thus far is applied by retaining the binary parameters used in the 
absence of a disk case. Accurate values regarding the minimum 
intervals f(emin) - To are given in Tables |6] and [T] for e - 0.01 d, 
while they are graphically demonstrated in Fig. |2] for the two 
more noise levels of e = 0.001 d and 0.03 d. 

Careful inspection of the results reveals that the O-C dif¬ 
ferences are detectable as long as the ML rate from L2 ex¬ 
ceeds -10“^' A/oyr^' for e - 0.001 d or -10“'° A/g yr“^ when 
the noise level takes higher values (i.e., e - 0.01 or 0.03 d). 
However, the O-C residuals produced by ML via the L3 point 
need a modestly longer baseline of observations compared to 
the corresponding process via the L2 point. The mass ratio q 
proves to be a factor that hardly influences the temporal range 
l(emin) - To, lying in a few decades no matter what the adopted 
noise level is. 


Table 6. Minimum time spans at e = 0.01 d when material, 
ejected from the more massive component, is carried away from 
the system through the L2 point (fully liberal case). 


^ = 0.5 

mL2 t- To 

[A/oyr-‘] [yr] 

' jlL 2 = 
9 

-lO-** A/oyr-‘ 
t —To 
[yr] 

- io -‘- 

2798 

0.9 

0.19 

26.8 

- 10 -“ 

885 

0.8 

0.18 

27.1 

- io -‘° 

280 

0.7 

0.18 

27.3 

- 10 -” 

88.5 

0.6 

0.17 

27.6 

- 10 -* 

28.0 

0.5 

0.17 

28.0 

- 10 -” 

8.8 

0.4 

0.16 

28.4 

- lo -*^ 

2.8 

0.3 

0.15 

28.9 

- 10-5 

0.9 

0.2 

0.13 

29.7 

- 10 -^ 

0.3 

0.1 

0.11 

31.1 


3.3. Case III: Non-conservative mass transfer 

We examine now what the combined effects of MT, ML, and 
AML processes might be in the modulation of an ETV diagram 
by distinguishing the case where the donor is the primary or the 
secondary component. Furthermore, we discern two subcases of 
each non-conservative MT case, addressing the general regime 
where a transient disk surrounds the gainer; the first concerns 
ML and AML realized by a hot-spot emission, while the second 
regards ML and AML from the L2/L3 point (non-captured mass 
from the gainer). 


3.3.1. Critical mass ratios when the donor is the more 
massive star 

If jS 6 [0,1] is the fraction of the MT rate gained by the accretor, 
then /3 - I represents the conservative case (where all the trans¬ 
ferred mass is captured from the gainer), while fi -0 represents 
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Table 7. Minimum time spans at e = 0.01 d when material, 
ejected from the less massive component, is carried away from 
the system through the L3 point (fully liberal case). 


Table 8. Critical mass ratios q„ as a function of the degree of 
liberalism fi when the donor is the more massive component (see 
text for details). 


5 = 0.5 

niLi t — To 

[Moyr-‘] [yr] 

Wls = 

q 

-lO-^Moyr-* 
r\ t-To 

[yr] 

-10-*' 

4173 

0.9 

0.20 

31.5 

-10-'* 

1320 

0.8 

0.21 

33.2 

-io-“> 

417 

0.7 

0.22 

35.4 

-10-® 

132 

0.6 

0.24 

38.2 

-10-** 

41.7 

0.5 

0.25 

41.7 

-10-’ 

13.2 

0.4 

0.28 

46.6 

-10-" 

4.2 

0.3 

0.31 

53.8 

-10-" 

1.3 

0.2 

0.36 

65.7 

-10-'* 

0.4 

0.1 

0.45 

92.0 


the fully liberal case (where all the transferred mass is rejected 
from the gainer). 

As far as the donor is the more massive component and pro¬ 
vided that Ml - -M 2 + 'h is valid, the MT rate captured by the 
secondary, M 2 , and the remaining fraction that evades the sys¬ 
tem, m, can be written as 

M2=-PMu (29) 


P 

ML through a 

^cr 

(r„rL = 0) 

hot spot 

^cr 

(r„ Tl t 0) 

ML through the L2 point 

^cr ^cr 

(rr, rL = 0) (r„ fL t 0) 

0.0 

0.72 

0.49 


0.1 

0.74 

0.53 


0.2 

0.77 

0.57 


0.3 

0.79 

0.62 


0.4 

0.82 

0.67 


0.5 

0.85 

0.73 


0.6 

0.88 

0.80 


0.7 

0.91 

0.88 


0.8 

0.94 

0.97 


0.9 

0.97 



1.0 

1.00 


1.00 


for q'cr are then presented in Table [8] Binary systems with mass 
ratio q greater than the listed values will reveal an increasing 
period (and a convex ETV diagram), as long as the driving pro¬ 
cesses are those examined in the current section. New q^^ values 
can emerge by substituting r 2 for r, in the permanent disk for¬ 
mation mode. 


and 


m = (l-yS)Mi. (30) 

In the first subcase (hot-spot-driven ML and AML from the 
secondary), we combine Eqs. (l29l l and (l30l l to express M\ and 
OThs as functions of the mass accretion rate M 2 - Then, Eq. O 
reduces into the following form: 


j = [-3(1 -q) - 3 yjqiq + l)ri- + | ' (? + 1 )'‘l 
1-/3 'iq^ + 2q-3 M 2 


(31) 


Proceeding to the second subcase (ML and AML from L2), 
the rates Mi and otl 2 are expressed as a function of M 2 as well, 
and finally Eq. ([TJ takes the following form: 


^ = [-3(1 -q) - 2> yjqiq + l)r, + | ' (9' + ^>1 

1 -13 3;L2(g+ l)^-3g^-2g ^M2 
P ' q+l ^M 2 ^ 

It is worth noting that on the r.h.s. of Eqs. OTI) and (l32l i. the 
expression within the brackets is a function of the mass ratio q 
alone. As a consequence, in the framework of our hypothesis set 
(and as long as p is treated as a free parameter), the mass ratio 
q becomes a crucial factor that modulates the variations (actu¬ 
ally the monotonicity) of the period function P(f), as well as the 
geometrical form (actually the curvature) of the respective ETV 
diagram. The critical mass ratio values q^r can therefore arise 
after determining the roots of this expression, for any /?, consid¬ 
ering the direct MT mode either in the presence of a transient 
disk orbiting the secondary or in the absence of a disk. This can 
be achieved by calculating r, from Eq. ([Thl l and tl from Eq. ([TST i 
in the former case and by setting = tl = 0 in the latter. Results 


3.3.2. Critical mass ratios when the donor is the less 
massive star 

As far as the donor is the less massive component and provided 
that M 2 = -Ml -H m holds, the following relations are valid for 
the mass accretion and the ML rates. Mi and m, respectively: 

Ml = -PM 2 , (33) 

and 

m = (l-yS)M2- (34) 

In the first subcase (hot-spot-driven ML and AML from the 
primary), we combine Eqs. (l33T l and (l34l i to express M 2 and mhs 
as a function of the mass accretion rate Mi and, subsequently, 
Eq.@ is written as 

^ = [-3(q' - 1) - 3 yj{q+ l)r, + | ' (9' + 1)'‘l 
P q+l ' M2 

In the second subcase (ML and AML from L3), both the 
accretion and the ML rate, M 2 and mL 3 respectively, are written 
as a function of Mi, modifying Eq. (|2]l as 


j = [-3(5' - 1) - 3 ^J{q+ \)r, + I ■ (9' + ^>1 

_ 3;L3(g+l)"-25-3 Ml 

P ' q+l M2' 


(36) 


The numerical expression within the brackets on the r.h.s. of 
Eqs. (I 35 I 1 and (l36l l depends exclusively on the mass ratio q of 
the system. Critical mass ratio values q^ then emerge by deter¬ 
mining the roots of this expression for different p values, con¬ 
sidering a transient disk around the primary, with r, and be¬ 
ing estimated via Eqs. ([TTI i and ([T9l l. respectively. Accordingly, 
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Table 9. Critical mass ratios as a function of the degree of 
liberalism (3 when the donor is the less massive component (see 
text for details). 


-210 -140 


t-T„ [yr] 

-70 0 70 


140 


210 


ML through a hot spot 
[3 qcT qcT 

(l, rL = 0) (r„ rL t 0) 


ML through the L3 point 

^cr ^cr 

(ri.,rL = 0) (rr,rLi*0) 


0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 



< 0.05 
0.05 
0.10 
0.17 
0.24 
0.32 
0.42 
0.53 
0.66 
0.81 
1.00 


<0.05 

<0.05 

0.07 

0.11 

0.17 

0.23 

0.31 

0.40 

0.52 

0.66 

0.83 


the qci values in the absence of a disk were derived by posing 
fr = fL = 0. Results for q^- values are then given in Table |9] 
Systems with mass ratios q greater than the listed q^r values 
are expected to show a decreasing orbital period (and a concave 
ETV diagram). The same procedure can be applied by substitut¬ 
ing r\ for Tr to find the corresponding q'cr values in the presence 
of a permanent disk. 

In Fig. [3] synthetic ETV diagrams have been constructed 
by means of Eqs. (fTSl l. dTSl) . and (l3^ for a q - 0.4 semi¬ 
detached system with a primary gainer of ri = 0.13 when the 
mass accretion rate is hxed at Mi = 10“* Mq yr"'. The geome¬ 
try of the anticipated ETV diagram for the same system is mod¬ 
ulated according to the MT mode and the degree of the liber¬ 
alism (i.e., the fi value). Indeed, for the system RR Dra, which 
is suspected of accommodating a transient disk orbiting the pri¬ 
mary star (Kaitchuck et al. 119851 1. a conservative MT rate close 
to 6 X 10"^ Mq yr"' is derived, while a value of 3 x 10"’ Mq yr"' 
is (under)estimated by the conventional approach (Zasche et al. 
12008b . Moreover, the observed period variations of X Tri could 
be ascribed to a non-conservative MT process with ML rates 
through L3 ranging from -2 x 10"^ to -6 x 10"*Moyr"' (see 
also Nanouris et al. 12013b . whereas the conservative approach 
fails to provide any reasonable explanation (Rovithis-Livaniou 
et al. [20b0b . 



Fig. 3. Synthetic ETV diagrams for aq - 0.4 semi-detached sys¬ 
tem when the donor is the secondary star (step s 1000 cycles). 
Squares and circles account for a conservative era {ji - 1.0) ei¬ 
ther in the absence or in the presence of a transient disk, respec¬ 
tively. Upwardly pointing triangles (J5 - 0.9) account for a hot- 
spot-driven non-conservative era in the presence of a transient 
disk, while crosses (fi - 0.7) and downwardly pointing triangles 
(fi - 0.5) represent a non-conservative era where material es¬ 
capes via the L3 point. (This hgure is available in color in the 
electronic form.) 

AML processes, the following differential equation is set and is 
solved for the orbital period P(t)\ 

P = (38) 

where the constant coefficient hg is given by the relation 
M 1 M 2 




5c5 


(Ml -H M2)'/3 


> 0 . 


(39) 


The solution of Eq. (l38T l under the initial condition (f. Pit)) - 
(To,Pfi is 


3.4. Case IV: Angular momentum loss through gravitational 
radiation 

The last studied individual case concerns the effects of gravi¬ 
tational wave radiation (GWR), a physical mechanism that is 
highly suspected to drive binaries evolution with mainly degen¬ 
erate components (e.g., Kramer et al. 120061 Antoniadis et al. 
12013b . In non-degenerate binaries, GWR seems to become a sig- 
nihcant source of AML once binaries reach the contact phase 
(Webbink [T976b . 

The rate at which a binary loses angular momentum because 
of the GWR mechanism is given by the well-known relation of 
Landau & Lifshitz (11962b : 


4 = - 


32 (InG) 

5^ 


7/3 


(Ml M 2 )’ 


(Ml -H M^WPP^ ’ 


(37) 


R(f) = 


Pi - \b^(t - To) 


(40) 


Then, integration of Eq. (|4]i under the initial condition (e, t) - 
(0, To) yields 




1 

f T 5 ) 

8 - 
5 





(41) 


Hence, through Eqs. (l40b and (l4Tb . the orbital period P(e) as a 
function of the continuous time variable e is 


(42) 


Pie) = - -hge 


Integration of Eq. (|5]l under the initial condition (e, AT(e)) 
(0,0) eventually gives the AT(e) function 


where c is the speed of light in a vacuum. By ignoring MT and AT(e) - ^ 
ML terms in Eqs. ([T]i and ^ and disregarding any other possible 8hg 




(43) 
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Table 10. Minimum time spans at e = 0.001 d and 0.01 d for sev¬ 
eral systems that lose angular momentum through gravitational 
radiation. 


t-To 

lyr] 

System 

Te 

[4] 

+ AI 2 

[Mq] e 

O-C noise level 
= 0.001 d e = 0.01 d 


1.00 

1 -tl 

933 

2951 


1.00 

5+5 

244 

772 


1.00 

8 -t8 

165 

521 

CC Com 

0.22 

0.74-0.4“ 

212 

670 

WUMa 

0.33 

1.24-0.6“ 

246 

779 

AH Cep 

1.77 

164-13* 

216 

684 

WR20a 

3.69 

83-h82“ 

135 

426 


°Hilditch etal. <19881 


^Burkholder et al. ( 119971 
"^Bonanos et al. ( 120041 

which practically pictures the theoretically anticipated form of 
the ETV diagram under the action of GWR-AML, and is proved 
to be concave for any cycle e, a result referring to the con¬ 
tinuously decreasing orbital period function Pie), expressed by 
Eq. dUi. 

The traceability of a 1 -i- 1 Mq, a 5 -(- 5 Mq, and a 8 -i- 8 Mq 
system are then examined, all with Pe = 1 d. As seen in Table [TOl 
even for a quite massive pair of 8 -(- 8 Mq, it is almost impossible 
to observe any GWR-driven O-C variations, since at least 160 
and 500 yr of dense monitoring are needed below noise levels 
equal to 0.001 and 0.01 d, respectively. In Table fTOl we also tab¬ 
ulate the corresponding time intervals f(emm) - Tq expected for 
the short-period contact binaries CC Com and W UMa (Pg < 
0.4 d) and, beyond these, for the two massive binaries AH Cep 
and WR 20a (Mi -i- M 2 > 30 Mq). Even for WR 20a, which is 
among the most massive binaries known to date (Bonanos et al. 
120041 1. the GWR effects on its ETV diagram would only become 
detectable after 135 yr of accurate recorded times of minima. 

An additional investigation is also performed to hnd out 
whether GWR-driven modulations are visible in a wide range 
of periods, varying from 0.1 to l.Od. All three hypothetical 
systems are explored, keeping in mind that a period threshold 
should exist for each case as a response to the contact conhg- 
uration (Eig. Si. The critical limit of the orbital period P^r was 
determined by considering radii of 3.6 P© and 5.2 P© for a 5 M© 
and a 8 M© dwarf star, respectively. Especially in the case of the 
1 -1-1 M© system, the critical period P^t was also determined con¬ 
sidering that one member of the pair is a white dwarf with radius 
equal to 0.01 P© (MSh-WD). 

A close inspection of Eig. |4] reveals that, as in the mag¬ 
netic braking case (see NKARLll), the shorter the orbital pe¬ 
riod, the shorter the time required. Close to the period thresh¬ 
old of each examined system, slightly less than two centuries of 
observations are needed for producing measurable O-C varia¬ 
tions. However, the MS-i-WD system seems to be an exceptional 
case that is likely to be highly detectable, demanding only hve 
decades of monitoring at most. This hnding mainly accommo¬ 
dates the group of eclipsing cataclysmic binaries with orbital pe¬ 
riods of a few hours (see, e.g., Pilarcik et al. 120121 and references 
therein). 

4. Discussion 

We examined the orbital evolution and the associated synthetic 
ETV diagrams of binary systems that experience some of the 


0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 4. Minimum time spans considering that gravitational radia¬ 
tion is the driving evolutionary process of a 1 - 1 -1 M© (squares), a 
5 -(- 5 M© (circles), and a 8 - 1 - 8 M© system (up triangles) with 
various periods for noise levels equal to 0.001 d (left vertical 
axis) and 0.01 d (right vertical axis). The labels depict the period 
thresholds at which each system attains the contact configura¬ 
tion. (This figure is available in color in the electronic form.) 


most common manifestations of MT mechanisms in short, hu¬ 
man life timescales. A complementary study was also carried out 
for GWR-induced modulations. The minimum range, required 
for the detectability of the ensuing O-C variations, was then es¬ 
timated following the methodology introduced in Nanouris et al. 
(1201 111 for a variety of both physical and observational parame¬ 
ters where the procedure had been proved to be sensitive. 

The conservative MT process was thoroughly investigated, 
and it was shown that for the direct-impact mode in the presence 
of a transient disk - which extracts large amounts of the orbital 
angular momentum - around the more massive component there 
is a critical mass ratio {q^ ~ 0.83) above which the period turns 
out to decrease, contrary to those expectations that emerged by 
the conventional approach. 

Different non-conservative MT modes were also examined 
by employing a parameter (3 to denote the degree of conser¬ 
vatism. Values of (3 close to zero account for violent non¬ 
conservative eras, while values close to unity account for highly 
conservative ones. TakingyS as a free parameter in their hot-spot- 
driven evolutionary models, van Rensbergen et al. (120081 1 have 
shown that short non-conservative eras soon after an RLOE ig¬ 
nition may take place with a variable (3, even reaching the value 
of (3 ~ 0.2. The present study revealed that in a non-conservative 
MT process induced by a hot-spot emission, AML is too weak 
to drive the orbital evolution, still allowing an increasing pe¬ 
riod even in strong but reasonable non-conservative eras (i.e., 
(3 < 0.4), however only for a narrow range of high mass ratios 
(i.e., q > 0.7) in systems where the donor is the primary star. 

Instead, ML via the outer points L2 or L3 seems to be a 
strong AML process, capable of reducing the period at a very 
rapid pace. Sytov et al. (I2007t have shown that the efficiency of 
accumulation in accretion disks of close binaries is not higher 
than 20-30% for typical values of turbulent viscosity, while it 
hardly reaches a rate of 50%, provided that the magnetic field is 
taken into consideration (Bisikalo 1201 Oil , impelling ML through 
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L3 at significant amounts. We showed that the period is expected 
to decrease even during weak non-conservative eras through the 
L2 point, while when L3 is the libration point, it was revealed 
that there is always a critical mass ratio above which the 
period is expected to decrease. The critical values decrease 
as the degree of liberalism is progressively widened. Unlike the 
hot-spot emission, this allows a period reduction even for very 
low mass ratios (i.e., q < 0.3) in systems where the donor is the 
secondary star. 

The above findings prove that ML through the L3 point 
is a forceful candidate mechanism for interpreting the concave 
ETV diagrams that represent semi-detached systems with the 
less massive star to be the contact component (see, e.g., Yang 
& Wei 12009] for an updated catalog of binaries that display this 
peculiar behavior). Among them, AF Gem, TU Her, AT Peg, TY 
Peg, Z Per, Y Psc, and X Tri are systems with low mass ratios, 
and consequently, they are strongly suspected to undergo a non¬ 
conservative MT era, accompanied by ML through the L3 point 
at a significant rate. Recently, light curves of several early-type 
systems with an accretion disk have been analyzed using the new 
improved code of Djurasevic et al. (12005112008II for estimating 
the physical disk parameters. Their analysis has indicated the 
presence of a bright spot, a region on the disk edge opposite 
(and close) to the L3 point, through which matter is expected to 
escape from the system. DL Cyg, p Lyr, AU Mon, V393 Sco, 
RY Set, and DQ Vel are low mass ratio systems that could expe¬ 
rience a non-conservative MT era. 

It is noteworthy that the MT process is capable of driving the 
orbital evolution in the limited temporal range of ETV diagrams 
at rates similar to the wind-driven ML counterparts (NKARLl 1), 
implying that these two mechanisms might be strongly com¬ 
petitive (see also Erdem et al. 120051 Erdem et al. 12007al bf. 
Consequently, ML through stellar winds should not be neglected 
in systems of semi-detached or contact configuration. Besides, 
the inferred critical mass ratios from the present analysis are ex¬ 
pected to change considerably when this kind of ML is involved 
in our calculations. Evidently, the simplistic p-q schemes seem 
to be insufficient in describing the short orbital evolution of a 
binary in which MT is not the leading evolutionary mechanism. 
This may be consistent with the findings of Qian (120011 1. who 
has shown that a critical value close to q^ ~ 0.40 appears to in¬ 
vert the monotonicity of the observed orbital evolution in contact 
binaries. However, this did not emerge clearly from the present 
study. The determination of the period monotonicity in the pres¬ 
ence of numerous physical processes of a different nature, al¬ 
though feasible, is rather complicated and requires a more care¬ 
ful investigation. 

The sufficiency of the commonly adopted parabolic repre¬ 
sentation of the long-term variations in ETV diagrams will be 
examined in a forthcoming paper through Taylorian expansions 
of the analytic AT(e) expressions. 
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